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ABSTRACT 

Stimulated  by  recent  work  by  Gregory  and  Qu,  it  is  shown  that  the  limit  of  local 
corner  cutting  is  a  continuously  differentiable  curve  in  case  the  comers  of  the  iterates 
become  increasingly  flat. 
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Local  corner  cutting  and  the  smoothness  of  the  limiting  curve 
/  Carl  de  Boor 

It  was  proved  in-ffify  that  corner  cutting  of  any  kind  converges  to  a  Lipschitz-continuous 
curve,  but  the  question  of  how  one  might  guarantee  that  the  limiting  curve  be  smoother 
than  that  was  not  considered  there.  Recently,  Gregory  and  Qu  [GQ]  took  up  this  ques¬ 
tion  and  established  sufficient  conditions  for  a  certain  systematic  and  local  corner  cutting 
scheme  to  give  a  limiting  curve  in  Since  [GQ]  use  the  same  parametrization  of  the 
successive  broken  lines  that  made  the  argument  in  [B]  so  simple,  I  became  intrigued  and 
took  a  look  at  what  one  might  say  in  greater  generality.  Specifically,  I  looked  for  conditions 
under  which  continuous  differentiability  of  the  limiting  curve  could  be  inferred  from  the 
fact  that  the  comers  of  the  broken  lines  flatten  out  eventually. 

It  is  the  purpose  of  this  note  to  prove  that  the  limit  of  any  ‘local’  corner  cutting 
scheme  is  in  C*Jlprovided  the  comers  of  the  broken  lines  become  increasingly  flatter.  A 
simple  example  is  given  to  show  that  this  condition  is  not  necessary,  while  another  simple 
example  shows  that,  without  ‘localness’,  the  condition  is  not  sufficient,  in  general.  Finally, 
as  an  application,  Gregory  and  Qu’s  nice  argument  in  [GQ]  is  redone.  k.p-  — 


1.  Cutting  corners 


In  this  section,  we  recall  the  setup  of  [B]. 

We  deal  with  a  sequence  (6„)£L0  of  broken  lines  in  which,  for  n  >  0,  bn  is  obtained 
from  bn- 1  by  a  ‘cut’,  i.e.,  by  replacing  a  curve  segment  by  the  subtended  secant  to  the 
curve.  This  means  that  all  the  vertices  of  bn  lie  on  bn- i,  i.e.,  bn  can  be  thought  of  having 
been  obtained  from  6n_!  by  interpolation.  This  observation  is  used  in  [B]  to  prove  that,  no 
matter  just  how  the  cutting  was  done  to  generate  the  sequence  (bn)  from  an  initial  broken 
line  &o  with  finitely  many  vertices,  b^  :=  limn_oo  bn  exists  as  a  Lipschitz-continuous  curve 
which  is  approached  uniformly  by  f>„,  i.e.,  lim„— oo  dist(6„,  boo)  =  0. 

The  argument  in  [B]  was  based  on  parametrizing  the  curves  appropriately.  If  (u,)  is 
the  sequence  of  vertices  of  bn  and  (#,)  is  a  corresponding  arbitrary  increasing  sequence  of 
numbers,  then  bn  can  be  parametrized  by 

(1.1)  bn(t )  :=  +  v*  v**  ,  <«-i  <<<<«,  all  i. 

Since  b„  is  obtained  from  bn-i  by  interpolation,  it  is  natural  to  choose  the  sequence  (t{) 
in  dependence  on  the  parametrization  of  6„_ j,  i.e.,  so  that  &„(/<)  =  v«  =  6n-i(*»)  for  all  t. 
With  this, 

bn  =  P n^n—li 

where  Pn  is  broken  line  interpolation  at  the  points  (<*).  Therefore,  ultimately,  bn  = 
P„  ■  •  • Pib0 ,  with  Pn-  ■■  Pi  a  linear  map.  Hence,  although  the  process  of  generating  the 
sequence  (&„)  is  nonlinear  (in  that  it  is  quite  arbitrary),  once  we  have  decided  on  how  to 
cut,  we  can  think  of  each  bn  as  a  linear  function  of  bo-  In  particular,  writing  &o  in  any  one 
of  many  reasonable  ways  as  a  sum 


of  scalar-valued  functions  ipi  with  vector  coefficients  Wi  €  IRd,  we  have 

bn  =  ^WiPn-'-Pitpi, 
i 

and  questions  of  convergence  or  of  smoothness  of  the  limit  can  be  settled  by  settling 
them  for  the  (presumably  simpler)  sequences  (Pn  •  •  •  Pi<p),  with  any  one  of  the  y?*.  A 
particularly  simple  choice  for  the  ifi  axe  the  truncated  powers  (•  —  Tj)+  (in  addition  to  the 
constant  function),  and  this  leads  to  the  conclusion  that  the  nature  of  the  limiting  curve 
can  be  understood  if  one  understands  what  the  particular  comer  cutting  process  does  to 
the  standard  corner,  i.e.,  the  broken  line  with  vertices  (0,0),  (1,0),  (2, 1). 

2.  Examples 


We  are  now  ready  to  consider  the  smoothness  of  the  limiting  curve,  having  understood 
that  it  is  sufficient  to  consider  the  case  that  bo  is  a  piecewise  linear  (real- valued)  function 
on  some  interval.  Then  each  bn  is  of  the  same  nature,  and  its  derivative,  dn  :=  b'n)  is  a 
step  function.  In  the  notation  adopted  in  the  preceding  section, 

(2.1)  dn  =  J^(-  -  ti)\  jumpj.  dn, 

i 

with  jumped  :=  d(t- f-)  —  d(t-)  the  difference  between  the  limit  from  the  right  and  the  limit 
from  the  left  at  t.  We  can  take  the  absolutely  largest  jump,  i.e.,  the  number  ||jumpQdn||oo, 
as  a  measure  of  the  extent  to  which  b„  fails  to  be  in  C1. 

ra 


(2.2)  Figure  The  change  in  the  derivative  as  the  result  of  a  corner  cut.  The 
two  areas  are  of  equal  size. 

It  will  be  useful  to  visualize  the  process  by  which  dn  is  obtained  from  dn-\.  Suppose 
that  bn  is  obtained  from  6„_i  by  replacing  f»n_i  on  [s,f]  by  the  linear  interpolant  to  6„_i 
at  s  and  t.  Then 

(d„  -  dn-i)(x)dx  =  (bn  -  bn-i)(t)  -  ( bn  -  6„_i)(s)  =  0. 
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Also,  dn  is  a  constant  (viz.,  the  difference  quotient  [s,  t]i„-i )  on  [a,#].  If  now  &„-i  has 
just  one  vertex  in  [a,t],  then  dn-i  has  just  two  steps  there,  hence  d„  —  d„_i  has  just  two 
steps  there  and,  as  f*(dn  —  dn_x)  =  0,  the  two  rectangles  which  make  up  this  integral  must 
balance;  see  (2.2)Figure. 

More  precisely,  let  J  be  the  jump  in  dn- 1  at  that  sole  vertex  in  [s,  t],  let  l  and  r  be  the 
parametric  distances  of  the  corner  from  its  left  and  right  neighbor,  and  assume  that  the 
two  new  vertices  (which  replace  the  vertex  being  cut  off)  occur  at  parametric  distances  //? 
and  ra,  respectively  (with  a, ft  G  [0, 1]).  Then  the  two  new  jumps  are  of  size  hiJ  and  hrJ, 
with  hi  +  hr  =  1  and  l/3hi  =.  rahr,  or, 


(2.3) 


hi 


ra 

l  ft  4-  ra 


1  -  hr. 


(2.4)  Figure  These  uniformly  nonsmooth  broken  lines  converge  to  a  smooth  limit. 

As  a  first  example,  we  show  that  the  maximum  jump  need  not  go  to  zero  for  the  limit 
function  to  be  C1.  As  shown  in  (2.4)Figure,  the  limiting  function  is  smooth  (it  is  the 
zero  function),  while  the  maximum  jump  in  the  first  derivative  stays  above  1  (in  absolute 
value).  Note  that,  in  this  example,  we  have  done  two  cuts  simultaneously,  i.e.,  in  terms  of 
the  setup  adopted  earlier,  we  are  showing  only  every  other  iterate.  Note  also  that,  strictly 
speaking,  each  of  our  cuts  involves  two  corners.  It  will  be  important  later  on  to  know 
that  such  an  example  (of  a  limiting  curve  being  smooth  even  though  the  absolutely  largest 
jump  in  the  derivative  of  the  iterates  is  bounded  away  from  zero)  can  also  be  given  when 
each  cut  involves  only  exactly  one  comer.  Such  an  example  can  be  supplied  by  applying, 
e.g.,  Chaikin’s  algorithm  to  the  initial  broken  line  in  the  above  example,  except  that  the 
points  on  the  two  segments  flanking  the  middle  segment  are  chosen  closer  and  closer  to 
the  farther  endpoint. 

Finally,  we  illustrate  the  fact  that  having  the  maximum  jump  in  the  first  derivative 
go  to  zero  is,  in  general,  no  guarantee  that  the  limiting  curve  is  C1.  For  that,  we  consider 
the  1-parameter  family  of  functions  fa  (see  (2.6)Figure)  given  by  the  rule  that 

(  0,  t  <  0; 

(2.5)  fa(t)  :=  |  at2,  0  <t<  ta] 

1 1  +  (<  -  l)2ata,  ta  <  t, 

with  ta  :=  1  —  y/l  —  1/a.  One  verifies  that,  for  o  >  1,  /<*  is  in  C1,  and  that  fa 
converges,  on  [-1,1]  say,  monotonely  and  uniformly  to  the  ‘standard  comer’  which  is  not 
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Cl.  The  actual  example  is  provided  by  replacing  selected  fa  by  interpolating  broken  lines 
of  sufficiently  fine  spacing  to  make  the  jumps  in  the  first  derivatives  as  small  as  one  pleases. 


(2.6)  Figure  These  smooth  broken  lines  converge  monotonely  to  a  nonsmooth  one. 
(The  circles  are  the  points  ( ta ,  fa{ta ))  of  (2.5).) 

Nevertheless,  if  the  comer  cutting  is  local,  then  having  the  absolutely  largest  jump  in 
the  first  derivative  go  to  zero  does  imply  that  the  limiting  function  is  in  C1.  This  is  the 
content  of  the  next  section. 


3.  Local  corner  cutting 

We  say  that  the  comer  cutting  is  local  in  case  any  cut  involves  exactly  one  comer. 
This  means  that  the  cut  endpoints  must  lie  in  the  interior  of  the  two  segments  which  form 
the  comer  being  cut.  Schemes  that  cut  all  comers  simultaneously  fall  into  this  category  as 
long  as  the  cuts  of  neighboring  comers  do  not  share  an  endpoint.  For  we  can  then  think 
of  them  as  having  been  carried  out  one  cut  at  a  time.  In  particular,  the  comer  cutting 
scheme  considered  in  [GQ]  is  local  in  this  sense,  as  are  the  schemes  considered  in  [R]  so 
many  years  ago.  It  follows  that  every  segment,  of  the  original  broken  line  as  well  as  of  any 
subsequently  generated  broken  line,  is  tangent  to  the  limiting  curve,  hence  the  situation 
depicted  in  (2.6)Figure  could  not  have  been  generated  by  local  comer  cutting. 

(3.1)Theorem.  The  limiting  curve  produced  by  a  local  comer  cutting  scheme  is  C1  in 
case  the  maximum  jump  in  the  Erst  derivative  of  the  iterates  b„  goes  to  zero  as  n  — *•  oo. 
The  converse  holds  in  case  b0  is  convex  and  for  arbitrary  comer  cutting. 

Proof.  It  is  sufficient  to  consider  the  special  case  that  60  is  the  ‘standard  corner’. 
Then  6o  is  convex,  hence  so  are  all  the  iterates,  with  6n  growing  uniformly,  and  pointwise 
monotonely,  toward  the  limit  function  &oo,  which  is  also  convex.  Let  b  be  one  of  the  iterates. 
Then  d  :=  b'  is  a  monotone  increasing  step  function.  Recall  from  (2.1)  that 

d  =  iumpe>rf» 

j 
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with  0  <  t\  <  . . .  <  tm  <  2  its  breakpoints.  We  consider  also  the  two  step  functions  d+ 
and  d",  given  by  the  rule 

m 

d±  =  ~  tiTi)+  d \ 

j=i 

with  to  :=  0  and  (for  the  sake  of  neatness)  tm+ 1  :=  2.  Then  jjd+  —  <  2jjjumpQ<iJ|00, 

while  d~  <  d  <  d+  pointwise,  since  d  is  monotone  increasing. 


(3.2)  Figure  Corner  cutting  contracts  the  ‘envelope’  formed  around  the 
derivative  d  by  the  step  functions  d~  and  d+. 

Further,  if  b*  is  obtained  from  b  by  cutting  off  exactly  one  comer,  and  d*  is,  correspondingly, 
the  derivative  of  b*,  then  (see  (3.2)Figure) 

(3.3)  d~  <  (d*)“  <  (d*)+  <  d+. 

This  implies  that  the  first  derivative  of  all  subsequent  iterates  lies  between  d~  and  d+ 
Hence,  if  HjumpQdnHoo  — ♦  0  as  n  —*  oo,  then  (d„)  is  a  Cauchy  sequence  in  the  complete 
normed  linear  space  of  all  bounded  functions  on  [0, 2]  (with  the  max-norm).  Consequently, 
dn  =  6'n  converges  uniformly  to  some  bounded  function  doo-  Now  consider  the  modulus  of 
continuity  Woo  of  this  limiting  function.  For  any  h  >  0  and  any  n, 

u 00(h)  =  sup  (d0o(t)  —  doc(s))  <  sup  (d+(t) -d~(s))  =:Un(h). 

0 <t-a<h  0 <t-s<h 

Note  that  wn  is  a  nondecreasing  step  function,  with  0  <  ZZ7n(0+)  =  supJ(d„(tJ+2+)  — 
dn(tj-i-))  <  3||jumpQdn||oo  — *  0  as  n  — ►  00.  Consequently,  for  every  e  >  0,  we  can 
find  n  and  6  >  0  so  that,  for  all  h  <  6,  u„(h)  <  e.  This  proves  that  u>oo(0+)  =  0,  and  so 
establishes  that  d oo  is  continuous. 

We  now  know  that  b„  converges  uniformly  to  some  Lipschitz-continuous  function  boo, 
while  dn  :=  b'n  converges  uniformly  to  some  continuous  function  d< ».  It  is  a  standard  result 
that,  therefore,  ^00  —  dooi  i-e-i  the  limiting  function  600  is  in  C1. 

For  the  converse,  assume  that  the  sequence  (6„)  of  convex  broken  lines,  all  defined  on 
the  interval  [0,2]  say,  converges  uniformly  to  some  function  b.  Assume  further  that  b  is 
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continuously  differentiable  at  the  interior  point  p.  This  means  that,  for  some  modulus  of 
continuity  u  (i.e.,  some  positive  function  oj  on  (0,oo)  with  a>(0+)  =  0), 

[t, p]b  :=  ( b(t )  -  b(p))/(t  -p)  =  b'(p)  +  0(w(|t  -  p\ )). 

Now  consider  J  :=  jumpp6„  for  some  n,  and  let  e  :=  \\b  —  6n||oo-  Then,  for  any  small 
positive  h, 

0  <  J  <  [p  +  h,p]bn  -  [p,p-  h]fen  <  [p  +  h,p]b  -  [p,p  -  h]b  +  4 e/h  =  0(v(h))  +  4e//i. 

Since  u;(0+)  =  0  and  e  =  ||6  —  6n||  — »■  0  Ets  n.  — >oo,  this  implies  that  J  must  be  small  when 
n  is  large.  4k 

By  going  to  a  larger  envelope,  the  argument  can  be  applied  to  somewhat  more  general 
corner  cutting,  viz.  one  in  which  no  cut  reduces  the  number  of  comers.  Since  each  cut 
generates  two  comers,  this  means  that  any  cut  is  restricted  to  cut  away  no  more  than  two 
corners.  To  handle  this  more  general  situation,  one  would  change  the  definition  of  d*  to 

m 

d±  =  -  <;* 2)+  jump,.  d, 

j= i 

with  further  adjustments  at  the  endpoints. 

(3.4) Corollary.  The  limiting  curve  produced  by  a  comer  cutting  scheme  which  never 
cuts  across  more  than  two  comers  is  C1  in  case  the  maximum  jump  in  the  first  derivative 
of  the  iterates  b„  goes  to  zero  as  n  — ►  oo. 

4.  The  Gregory-Qu  result 

As  an  application  of  (3.1)Theorem,  we  now  consider  the  Gregory-Qu  scheme,  in  which 
bn  is  obtained  from  6„_j  by  a  simultaneous,  non-interfering,  cutting  of  all  corners;  hence 
the  scheme  is  local  in  the  sense  defined  earlier.  In  fact,  [GQ]  assumes  that  the  new  vertices 
generated  are  in  the  interior  of  old  segments  (i.e.,  that  all  a  and  /?  are  positive,  in  the 
notation  used  there  and  introduced  below).  But  if  we  allow  also  trivial  cuts,  i.e.,  cuts  that 
begin  and  end  at  the  same  vertex  (a  =  0  or  /?  =  0),  then  this  scheme  models  any  local 
corner  cutting. 

To  prove  that  the  limit  is  in  C1,  it  is  therefore  sufficient  to  prove  that  the  jumps  in 
the  first  derivative  go  to  zero.  For  this,  we  discuss  the  scheme  in  terms  of  the  step  function 
which  is  the  first  derivative  of  the  broken  line  in  question.  A  look  at  (4.1)Figure  might  be 
helpful. 

The  single  jump  of  height  J,  with  left  and  right  segments  of  length  l,  r,  spawns 
two  jumps,  a  left  one  of  height  J/ij,  with  hi  :=  rar/(l0i  +  rar)  and  with  segments  1/  := 
1(1  —  a/  —  0i)  and  r/  :=  //?/  +  rar,  and  a  right  one  of  height  Jhr,  with  hr  :=  1  —  hi  and 
with  segments  lr  :=  n  and  rr  :=  r(l  —  ar  —  0r).  Since 

h  1  1  _i  1 

ni  1  +  {0t/ar)l/r  r  1  +  (ar//?,)r/Z  ’ 
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(4.1)  Figure  The  change  in  the  derivative  due  to  one  step  of  the  Gregory-Qu 
process. 


we  can  be  assured  that  ht  and  hr  are  uniformly  smaller  than  1  (hence  the  limiting  curve 
is  in  C1)  provided  we  can  show  that  the  local  mesh  ratio  l/r  is  bounded  away  from  0  and 
oo.  For  this,  consider  the  local  mesh  ratios  UJn  and  lr/rr  spawned  by  the  cutting  of  this 
corner.  We  find 


l  l/r  i  = 


/(I  -  at-  fr) 

101  +  rar 


#), 


with 


m  :« 


t(  l-ai-ft) 

t0i  +  ar 


This  function  is  increasing  on  [0,  oo),  starting  at  0  with  a  value  of  0  and  a  slope  of  L'{ 0)  = 
(1  -  a,  -  fa)/ ar  at  0  and  taking  the  limiting  value  L(oo)  =  (1  -  a,  -  0t)/0l.  Consequently, 
L  maps  the  interval  (0,oo)  into  the  interval  [0,£(oo)].  Further,  if  L'{ 0)  >  1,  then  L  has 
an  attracting  fixed  point  in  that  interval,  viz.  the  point  F  :=  (1  —  a/  -  0i  -  ar)/0i  — 
(I'(0)  -  1  )ar/0i.  This  means  that  L  maps  any  interval  [ 1/M,M }  containing  F  (and 
contained  in  [0,  oo])  into  itself. 


(4.2)  Figure  The  function  L  contracts  around  the  point  F. 
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By  symmetry,  rr/lr  =  R(r/l),  with 

R(t)  := 


i(  1  -a  r-  0r) 
01  +  tar 


a  function  which  maps  [0,  oo)  into  the  interval  [0,  K(oo)],  and  which  has  the  fixed  point 
G  :=  (1  —  ar  —  0T  —  0{)/ar  =  (R'( 0)  —  1  )0i/ar  in  case  R'(0 )  =  (1  -  ar  —  0r)/0i  is  greater 
than  one.  This  means  that  R  maps  any  interval  [1/M,  M]  containing  G  (and  contained  in 
[0,  oo))  into  itself. 

We  conclude  that  the  local  meshratios  are  bounded  away  from  0  and  oo  provided  the 
fixed  points  F  and  G  are  eventually  bounded  away  from  0  and  oo.  As  [GQ]  point  out,  this 
can  be  guaranteed  by  having  a  and  0  eventually  bounded  away  from  zero,  i.e.  having  both 
a  :=  liminfa  and  0_  :=  liminf/3  be  positive,  and  having  2/(0),22'(0)  >  1  when  formed 
with  ai,ar  =  a  :=  limsupa  and  0i,0r  =  0  :=  limsup/3.  This  amounts  to  the  conditions 

0  <  Q.,0^  and  a,  0  <  1  -  a  —  0. 

When  these  conditions  are  violated,  we  cannot  be  certain  that  the  local  meshratios 
stay  away  from  0  or  oo,  hence  the  reduction  factors  hi  or  hr  may  come  close  to  1.  This 
does  not,  of  itself,  imply  that  the  limiting  curve  has  comers.  But  the  above  discussion  is 
sufficient  to  show  that  the  limiting  curve  has  corners  if  L'{  0)  or  21' (0)  are  uniformly  below 
1. 

We  discuss  this  only  for  the  case  of  constant  a  and  constant  0.  Assume,  for  example, 
that  a  and  0  are  such  that 

#(0)  - - <  I- 

Then,  starting  with  the  ‘standard  comer’,  the  cutting  process  generates  a  sequence  of 
vertices  proceeding  to  the  right  with  associated  local  mesh  ratios  r/l  equal  to 

R(l),R\l)  =  R(R(l)),R3(l),... 

which  decay  geometrically  to  zero.  In  fact,  22"(1)  ~  (i?'(0))n  asn  — *  oo.  The  corresponding 
reduction  factors  therefore  satisfy 

(4-3)  =  l+wimi) =  1 "  + ~ 1  -  («/«*'(»))” 

We  want  to  show  that  the  corresponding  sequence  of  jumps  is  bounded  away  from  zero. 
Since  this  is  a  decreasing  sequence,  it  is  sufficient  to  show  that  its  limit,  the  infinite  product 

fK">. 

n=l 

is  positive.  This  is  the  same  as  proving  that  the  infinite  series 


5>M”> 
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is  finite.  From  (4.3),  lnh^  ~  —(a/ ft)(R'( 0))",  i.e.,  the  terms  of  the  sum  behave  like  that 
of  a  convergent  geometric  series,  hence  the  series  converges. 


The  foregoing  analysis  also  explains  the  fractal  nature  of  the  resulting  curves  (when 
using  constant  a  and  ft).  For  it  shows  that  the  height  of  a  particular  jump  or  the  meshratio 
at  a  particular  breakpoint  of  the  nth  iterate  is  the  result  of  two  contending  fixed  point 
iterations,  L  and  R,  with  the  influence  of  each  entirely  determined  by  the  particular  se¬ 
quence  of  right  and  left  turns  taken  to  reach  the  breakpoint  in  question  from  the  original 
breakpoint.  In  particular,  we  expect  any  collection  of  jumps  sharing  the  first  few  of  these 
turns  to  look  like  any  other  collection  of  jumps  sharing  the  first  few  of  these  turns. 

It  would  be  interesting  to  explore  further  special  situations  when  the  fixed  points 
coincide.  For  example,  both  fixed  points  (for  L  and  i2)  are  1  exactly  when  a  +  ft  =  1/2. 
In  this  case,  all  the  mesh  ratios  are  the  same,  hence  the  left  factors  hi  are  all  the  same 
as  are  all  the  right  factors  hr.  The  requirement  that  hi  =  hr  is  satisfied  exactly  when 
a  =  ft.  Thus  the  imposition  of  both  requirements  leads  to  a  =  1/4  =  ft  which  is  Chaikin’s 
algorithm. 
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